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POTENTIALS 
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Saint-Petersburg (Leningrad) State University, 198504 St. -Petersburg, Russia 

Multidimensional factorization method is formulated in arbitrary curvilinear 
coordinates. Particular cases of polar and spherical coordinates are consid- 
ered and matrix potentials with separating variables are constructed. A new 
class of matrix potentials is obtained which reveals a double degeneracy or 
equidistant splitting of energy levels (hidden symmetry). 



1 Introduction. 

The Factorization Method is an effective tool in search and construction of exactly solv- 
able Schrodinger problems [1] and of the models with equivalent (almost coinciding) energy 
spectra [2], [3]. The mainstream of the method was related with one-dimensional problems 

Recently, the multidimensional generalization of the Factorization Method was elaborated 
[1] - [6], and its supersymmetric origin was demonstrated [7], |8]. The main results for 
multidimensional case were obtained in Cartesian coordinates. Meanwhile, the symmetry 
properties of original potential mark out the preferred coordinate systems related to orbits 
of the symmetry group (separation of variables). Usually, these coordinate systems are 
curvilinear ones. 

Thus, it would be interesting to reformulate the Multidimensional Factorization Method 
in curvilinear coordinates and to diagonalize matrix potentials appeared under the multi- 
dimensional Darboux transformation (Sect. 2). Thereby, starting from scalar potentials 
with separated variables a class of matrix potentials will be found which is also amenable 
to separation of variables . The polar and spherical coordinates and the particular case of 
Coulomb potential will be considered below in detail. 
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The full variety of Hamiltonians interrelated by the intertwining relations obeys twofold 
degeneracy of the spectrum. The latter has [2], [3], [7], [8] the supersymmetric interpretation: 
spectra of systems with even and odd fermionic modes coincide. 

However, the Factorization Method allows the construction of potentials with additional 
twofold level degeneracy without a simple supersymmetric interpretation. Together with 
standard supersymmetry of the full system of equivalent Hamiltonians, this means fourfold 
degeneracy of levels. Such potentials are built in Sect. 3, and the Factorization Method 
in curvilinear coordinates is more convenient for their description. More general class of 
potentials whose energy levels are paired with a given constant (equidistant) splitting is 
found. 



2 Factorization Method in curvilinear coordinates. 

We shall start from the case of two-dimensional space with arbitrary curvilinear coordinates 
[q^, q"^), where the interval ds"^ = gikdq^dq^ is given by the corresponding metric tensor Qik- 
The required generalization of Factorization Method includes: 

1) constructing the operators Qf {I = 1,2), which factorize the initial Hamiltonian: 

2 

2) constructiing the matrix Hamiltonian 

^?i" = 4^ + v;»>, (1) 

which is connected with H^^'^ by intertwining relations; 

3) constructing the scalar Hamiltonian 



i/(o) = _lA + \/(°), 
2 



also intertwined with H-f^ 



(1) 

ik ■ 

These necessary steps of the Factorization Method provide the well known links [1] - [6] 
between the spectra of operators H^^\ hII\ H^"^^ and between the corresponding eigenfunc- 
tions ^(1), 

The expression for the operators can be found by using the representation of the 
Laplace operator (kinetic terms in the Hamiltonians above) in arbitrary curvilinear coordi- 
nates: 
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where g'^'' is a contravariant tensor inverse to the metric tensor Qik, and the covariant deriva- 
tive operator Vj is defined by means of three-index Christoffel symbols F : 

d 

Viip = diip = 7--V? ('^ — scalar), 

ViAj = diAj — r^j^fc {Aj — covariant vector), 



ViA^ = SiA-'+r^.A" {A^ - contravariant vector), (2) 



ViBji = diBji — T'^iBki — ^uBjk {Bji — second rank covariant tensor). 



etc. The Christoffel symbols are expressed in terms of the metric tensor Ff^ = ^g'^y^p- + 

— ; and they are introduced in differential calculus of tensor objects to take into 

account the change of basic vectors and metrics under parallel transport . In a particular 
case of Cartesian orthogonal coordinates = 6ij, Tf- = 0, Vi = dj. 



Let us choose the operators Qf in the form: 



Qf.-L(TV, + (V«))=i^(TV, + (a«)). 

where exp(— x)— one of the two solutions of the Schrodinger equation H'^^^ exp (—x) 
Eexp (— x), and 



(3) 



Q'^ = 9''Qt, [gr> = (V.V.x), [Q7,Q,] = [Qi,Qt] = 0. 

The presence of the volume factor ^/det(gik) in the scalar product and the property Vigkj = 
are necessary to prove that operators Qf" are mutually conjugat^. 
Then, i/^o) can be written in a factorized form: 



(V,x)(V'x)-(V/V'x) 



+ E, (4) 



where the covariant derivative in Qj^ acts onto a scalar Vi\l/ = di"^, and covariant derivative 
in Q'l~ acts onto a vector in accordance with ([2]). 

To build the matrix Hamiltonian H^^^ let's use the operators Pj^ = ^/gtikQ^^ , P^^ = 
gi^p^ _ ^^ikQT (^g ^ det (yfjfc, with eik being completely antisymmetric unit pseudotensor) . 

The important property of orthogonality PkQ^^ = PkQ^^ = Qt^^^ = Qk^^^ = 
guarantees, in particular, intertwining the operator and the operator Hf'''^ : 

H^^^'' = QlQ^^^ + P'=+ + 5^E = + {ViV\) (5) 



^We remind that the flat space only is considered in this paper. 
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Similarly to the case of Cartesian coordinates |1] - [6] , these relations lead to the connection 
of spectra of iJ*^"-* and H^-^^ : {En^} C {En^} and of their wave functions: 



,Q,¥'\E), ¥'\E) 



\/E - E \/E-E 
Analogously, the operator h[^^ is intertwined with the scalar Hamiltonian: 

+ E = ij(°) + Ax, (6) 



2 2 



(V,x)(V'x) + (VzV'x) 



whose spectrum {En^} also lies in the spectrum of H^^^''. An arbitrary point of the spectrum 
of H'^^^ coincides either with En^ or with Ei^^ : 



=p^¥^\e)- ¥^\e) 



^pk + ^,W^E) 



\/E-E " ' ' VE-E 

For illustration we consider the case of polar coordinates (p, ip) on the plane: 



9ij 



( p2 ) ; - p; r 



2 _ p2 
12 ~ 21 



P 



From 



one obtains the explicit form of potentials {di = d/dp, 82 = d/dip) 



y(o) 



p'^ p"^ p 



+ E, 



(7) 



dAx'-Ax v(^^ + j,dlx + -Ax 



P 

{diX? + -,{d2X? + + + -{d,x) 

p'^ p^ p 



E. 



(9) 



For the particular case of centrally symmetrical x = x{p)i the matrix potential is diagonal 
just in polar coordinates: 



(i)fe 



V^^\p) + dlx + E _ 

V^^\p) + \d^x + E 



The method can be generalized to a space of arbitrary dimension. For the case of physi- 
cally interesting three-dimensional space, the operators are second rank tensors: 
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The intertwining relations are: 



(i)fc 



(2)fc 



(l)m 



o-(2)fcp+ 
(2)fc. 



p+ rr(l)m. 

r(2)fe. 



where 



/J(0) 



(V,x)(V'x)-Ax 



h\'^^ = Qrg'=+ + p-p'^+ + £; = 5^//(o) + (V,V'=x); 

(V;x)(V'x) + Ax 



^(3) 



2 2 

For the particular case of spherical coordinates (r, 9, (/?), one has: 

5-11 = 1; 6-22 = 5-33 = sin^ 6', 



rL = = -r sin^ ^; = - sin ^ cos e- vl^ = rfg 



r^3 = cote, 



and 



(aix)^ + -,{d,xf + -T^Mxf - dh - \x - -Ax - ^a.x - 



+ E, 



(l)fc 



sin^ 6* 



+ dlx ^ (^dAx - 
did2X--Ax V^^'^ + ^dlx + ldrX 
did,x--Ax ^{d^d^x-^oiOd^x^ v^(o) + _L^a|^ + iaiX+^a2X 



d2diX - cot ^^ax 



(2)fc 



(l)fc 



X — )■ — X 
y(o)(^)^y{3)(^) y ' 



From these expressions, one can conclude that in terms of the spherical coordinates in 3- 
dimensional space the matrix potentials V^^\ V^'^'> are also diagonal for spherically symmet- 
rical case X = x(^)- Thus, the scalar problem amenable to separation of variables produces 
the matrix problems which allow the separation of variables as well. 
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From this point of view, the problem of the Coulomb potential and its Darboux trans- 
formation, considered earlier in [9], looks interesting: 

f 0\ /+-00\ 

Vl^^^ = o' ; V/'^' = o' . (10) 
\000/ \000/ 

Thus, the Factorization Method developed above in curvilinear coordinates allows to un- 
ravel the structure of the matrix Coulomb potential. When interpreting as a vector 
particle wave function, then it follows from (ITU]) that bound states in matrix Coulomb po- 
tential exist for one polarization (along r) only, and the motion with two other polarizations 
is free. 



3 Matrix potentials with pairing of levels. 



The interrelation of spectra of the Hamiltonians if for arbitrary space dimension d was 
studied in [1] - [6], and it was proven there that the spectrum of the joint Hamiltonian 



H 



/ ij(o) 








\ 





iJ(l) 













... 




V 





H^'^^ 


/ 



(11) 



is twofold degenerate (may be, excluding the ground state). This degeneracy has [7], [8] 
a supersymmetric interpretation: the spectra of Hamiltonians coincide for even and odd 
number of fermionic excitations: 

/ i/W \ / if(i) 

Hb= \ ; Hf= i ii(3) 

\ ... J \ ... 

At the same time, a class of potentials with higher degree of degeneracy exists. In order 
to construct the potentials with additional degeneracy in 2-dimensional space, we examine 
the case with identical potentials V^^^ and V^'^^ : 5V = F^^) _ •^/(o) = Ax = 0. The real 
solutions of the Laplace equation are simply expressed in Cartesian coordinates in terms 
of arbitrary analytical functions of complex variable z = x + iy : x = F{z) + F{z) = 
2ReF{z), z = x — iy. 

However, polar coordinates are more convenient to study the physical properties of po- 
tentials V^"'\ The general solution of the Laplace equation, nonsingular for finite p, is: 



X 



^ amXniip,'^), 



Xr, 



p™ sm.{mip + 5r, 



(12) 



m=0 
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where am, Sm are arbitrary real numbers. Substituting these solutions into dZ]) - (j9]), one 
obtains a class of nonsingular potentials: 

_ oo oo 



2 

m=l n<m=l 



(i)fe 



where 



V/^^ ^ Z]m=2 "^("^ ~ l)ttmP'^ ^ COs(mV9 + 6m] 

[m - l)amp'^~^ cos(mv9 + 6m) V2^^ ^ 

(14) 

oo 

yW 1 _ y(o) ^ ^ _ i)a^p"^-2 sin(m¥; + 5^); 



m=2 

oo 



yii)2 ^ yio) _ J2 m{m - l)a^p™-2 sin(m(^ + 6m). 



m=2 



The asymptotic behavior of (fT3!) and (fT^ for large p (and arbitrary </)) is determined by 
the first term in (IT^ . and it guarantees the existence of bound states. The particular cases 
of potential V^'^^ of this kind are the harmonic potential V^'-^^ = ap'^ and the anharmonic 
potential of the form V^^^ = ap^ + /3p^ + Aa(5p^ cos(v3 + 6). 

Due to coincidence of potentials V^'^'^ and V^'^\ the wave functions of matrix Hamiltonian 

(1\ k 

HI'', built by means of 2-dimensional Darboux transformation from the wave functions of 
scalar Hamiltonians iJ'-^^ and H^'^\ read, 

v^a) = ^ Qi^^'\E(:^), = ^ Pr¥^\E^:^ ^ E(°)), 

and they have the same eig envalue E^r^^ = E^^ = E^\ In other words all energy levels in 
potential V^*^^'*'^ are twofold degenerate. Coincidence of these levels with levels of H^^'> and 
means that the joint Hamiltonian 

/ 

H=\ 

V 

has fourfold degeneracy of all levels. 

A more general problem corresponds to 6V = V^"^^ — 1^'-°^ = e = Const. In such a case, 
evidently En = En^ + e, leading to related equidistant splitting of all levels of the matrix 
Hamiltonian i/f^' : 

_ ^(0) ^(1) _ ^(2) _ ^(0) , 
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The explicit form of such potentials is obtained by solving the equation Ax = e : 

^ oo 

where Xm is given by f|T2l) . Then, 

1 1 °° 

8 2 ^-^ 

oo ^ oo ^ 

+ ^ mnamanP"'^""^cos[(m - n)v3 + 5m - + ^ mamp'"sin(mv? + (5m) + 

n<m=l m=l 

^(1) ^ ( ^1 ^ Z]m=2 "^("^ ^ l)ar„P"'~^ COs(mV9 + Sm] 

\ YZ.=2 ~ l)"mP"'"^ COs{mip + 6m) ^^2 ^ 

with 

oo 
m=l 

and 

oo 

Analogously, in three-dimensional space one can consider potentials with condition, 

SV = - = Ax = 0. (15) 
Its solutions are expressed in terms of spherical harmonics: 

oo 

j=0 

j 

Y^{e, y^) = J2 c^mPpicose) sin(m(^ + 5J, j = 0, 1, 2, ... 

These functions allow to reconstruct the set of potentials V^^\ ^/^'"^) ^'•'^^ with spe- 

cific properties of spectrum degeneracy. Usually, in the framework of Factorization Method 

/ i/W \ ( \ 

4x4 Hamiltonians Rb = ( g ^(2) 1 Hp =1 q ^(3) 1 are equivalent: their 

spectra coincide (supersymmetric degeneracy [7], [8]). The condition ( ITSj) leads to twofold 
degeneracy of a part of levels in Hb and therefore, in Hp, i.e. to fourfold degeneracy of a 
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part of levels in the joint Hamiltonian (ITT]) . Generally speaking, this additional degeneracy 
appears in three-dimensional case for a part of spectrum only, since some additional energy 
levels (common for Hf'^^ and Hf''^) exist, which have no direct relation to the condition 
5V = 0. 

Note added in translation. We are grateful to our colleagues and, especially, to Juan Mateos 
Guilarte, for their interest to this old paper and for encouragement to translate it into 
English. 
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